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A Prediction  Interval  Approach  to  Defining  Variables  Sampling  Plans  for  Finite 
Lots  Required  to  be  of  High  Quality:  Single  Sampling  for  Gaussian  Processes 

by 

Kenneth  K.  Fertig  and  Nancy  R.  Mann 
ABSTRACT 

A prediction-interval  approach  and  an  assumption  of  a Gaussian  manufacturing 
process  are  used  to  derive  a variables  sampling  plan  applicable  to  finite  lots 
required  to  be  of  very  high  quality.  Values  tabulated  for  calculating  accept- 
ance regions  are  such  that,  with  high  probability,  accepted  lots  will  have  zero 
defects.  Further,  tables  are  given  for  selecting  sample  size  for  a specified 
combination  of  lot  size  and  acceptable  quality  level.  Comparisons  show  that, 
for  a fixed  risk  level,  substantial  savings  in  required  sample  size  can  be 
effected  over  those  specified  by  comparable  hypergeometric  sampling  plans. 

Lots  ranging  in  size  from  5 to  100  are  considered. 

Bounds  for  the  probability  of  accepting  a lot  with  a fixed  number  of  defects 
and  coming  from  a prescribed  manufacturing  process  are  derived  and  prove  to 
be  very  tight.  Not  only  are  these  bounds  useful  in  defining  the  required 
sample  size  in  the  sampling  plan  context,  but  the  method  of  derivation  has 
application  to  other  areas  of  sampling,  e.g.  from  truncated  or  stratified 
distributions. 
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1 . INTRODUCTION 


In  the  following,  consideration  is  given  to  sampling  from  finite  lots  for 
which  the  underlying  process  distribution  is  Gaussian  with  unknown  mean 
and  variance.  It  is  assumed  that  accepted  lots  are  required  to  be  of 
very  high  quality,  that  is,  with  high  probability,  accepted  lots  should 
have  zero  defects.  An  item  will  be  said  to  be  defective  if  its  measured 
parameter,  x,  is  less  than  some  specified  value  L (greater  than  a specified 
U). 

Suppose  one  were  to  make  no  attempt  to  use  the  knowledge  concerning  the 
distributional  form  of  the  process  that  generated  the  lot.  Acceptance 
or  rejection  of  the  lot  would  then  be  based  solely  on  the  attributes  data, 
i.e.  whether  or  not  the  observed  x’s  were  less  than  L (or  exceeded  U) . A 
required  sample  size  would  probably  be  determined  for  a specified  combina- 
tion of  lot  size,  rejectable  quality  level  and  acceptable  quality  level 
by  use  of  the  hypergeometric  distribution,  perhaps  from  tables  of  Lieberman 
and  Owen  [8],  Use  of  binomial  sampling  plans  which  inherently  make  no 
provision  for  the  finite  size  of  the  lot,  can  lead  to  sample-size  require- 
ments that  are  unnecessarily  large,  perhaps  even  exceeding  the  lot  size. 

Here  it  is  assumed  that  we  record  the  variables  data,  the  sample  measure- 
ments, Xj,  x^,  of  the  normal  variate  that  determines  whether  or  not 
an  item  is  defective.  One  would  expect  that  the  required  sample  size 
associated  with  a specified  combination  of  lot  size  and  acceptable  quality 
level  based  on  the  additional  information  x^,  ...,  xn  and  an  assumption 
of  normality  for  the  process  distribution  would  be  less  than  required 
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solely  on  the  basis  of  whether  or  not  the  x's  are  less  than  I.  (greater 
than  U) . What  is  desired,  therefore,  is  a small-lot  variables  sampling 
plan  (one  that  makes  use  of  sample  measured  values  and  the  distributional 
form  of  the  process)  that  can  be  applied  to  lots  required  to  be  of  very 
high  quality.  In  the  following,  only  a lower  specification  limit  L will 
be  discussed.  The  theory  developed  can  be  applied  to  upper  specification 
limits  by  appropriate  changes  in  notation  and  sign. 

The  small-lot  Gaussian  sampling  plans  of  Fertig  and  Mann  [3],  the  theory 
for  which  is  described  in  [5],  are  based  on  a decision-theoretic  approach 
and  an  economic  loss  function.  Use  of  the  loss  function  requires  that 
rejected  lots  be  screened  and  defective  items  be  replaced  by  nondefectives. 
Because  of  this  requirement,  these  sampling  plans  cannot  be  applied,  for 
example,  to  situations  of  destructive  testing.  The  sampling  plans  to  be 
described  in  the  sequel,  however,  do  not  depend  on  the  utilization  of  a 
loss  function  and  can  be  used  even  if  the  sampling  situation  involves 
destructive  testing.  Moreover,  these  plans  are  based  on  what  will  generally 
be  considered  a more  classical  approach  involving  none  of  the  Bayesian 
methods  required  in  the  decision  theoretic  approach,  e.g.  the  use  of 
"optimal"  priors  on  the  process  parameters  to  compute  the  posterior 
expectation  of  the  number  of  defects  remaining  in  the  lot  after  sampling. 

The  new  sampling  plans  are  "Double  Zero"  sampling  plans  in  the  sense 
described  by  Ellis  [ 2] , that  is,  lots  are  required  to  have  zero  defects 
with  high  probability  and  detection  of  a defect  in  a sample  requires 
rejection  of  the  corresponding  lot.  Ellis  [2]  maintains  that  such  plans, 
which  are  used  currently  at  Pratt  and  Whitney  Aircraft,  "are  especially 
effective  in  the  metal  working  industry."  The  new  plans  differ  from  Ellis' 
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in  that  these  are  based  on  actual  variables  data  while  his  are  not  and 
here  there  is  no  consideration  given  to  Ellis’  ’’gray  areas”  of  marginal 
quality. 


-3- 


2.  BASIS  OF  THE  SAMPLING  PLANS 


A sampling  plan  can  be  described  in  terms  of  its  acceptance  criteria  and 
the  manner  of  specifying  the  sample  size  as  a function  of  lot  size  and 
given  requirements  on  quality  level.  Here,  the  ..pproach  used  in  defining 
the  acceptance  criteria  is  based  on  the  theory  of  classical  prediction 
intervals . Thus,  for  a sample  of  size  n from  a lot  of  size  N,  the  sample 
variables  data  x^,  ...,  x^,  are  used  to  predict  whether  or  not  the  N-n 
unsampled  items  in  the  lot  are  all  nondefective.  The  actual  prediction 
interval  is  computed  on  the  basis  of  the  sampling  distribution  of 


R ifX-Y,  )/S 
n,N-n  - l,N-n 


(2.1) 


where  Y,  ..  is  the  smallest  observable  value  of  the  random  variate  X in 
l,N-n 

the  unsampled  portion  of  the  lot  and  X and  S are  the  random  variates  with 

realizations  x and  s,  respectively,  calculated  from  Xj,...,  x^.  If  we  let 

8 be  the  lOOyth  percentile  of  R .,  , then  (X-Q  ,,  S,  + ") 

Y,n,N-n  p n,N-n  y,n,N-n 

is  the  lOOy  percent  lower  prediction  interval  for  Yj  N . 


Methodology  for  relatively  easy  and  inexpensive  computer  calculation 
of  9^  n ^ n was  developed  by  Fertig  and  Mann  (6).  Ihe  procedure 
depends  on  the  numerical  integration  of  a single  integral  (shown  in 
Section  5),  in  contrast  to  the  earlier  approach  of  llahn  [7]  which 
involved  a multivariate  t-distribution  and  numerical  integration  of  a 
double  integral. 


Care  must  be  exercised  in  the  interpretation  of  a prediction  interval. 

Specifically,  in  this  case  it  means  that  if  pairs  of  samples  of  sizes  n 

and  N-n  are  created  repeatedly  by  the  manufacturing  process,  the  100  Y 

percent  random  prediction  interval  (X-0  ..  S,  * ")  with  X and  S 

1 r y,n,N-n 


computed  from  the  first  sample  of  size  n of  the  pair,  will  contain  N n 
the  smallest  observable  variate  in  the  corresponding  second  sample  of  the 
pair,  100Y  percent  of  the  time.  If  the  prediction  interval  contains  the 
smallest  observation  y.  vl  , then  of  course,  it  contains  all  the  observa- 
tions  in  the  second  sample,  the  unsampled  portion  of  the  lot  in  the 
present  context. 

2.1  Acceptance  Criteria 

For  the  sampling-plan  tables  presented  in  this  paper,  we  have  set 

Y s 0.D0.  The  value  of  9 VI  for  each  n and  N considered  is  given  in 
1 Y>n»N-n 

Table  1.  The  basis  for  choosing  the  sample  size  n will  be  discussed  in 
section  2.2.  Once  the  sample  ...»  xR  has  been  observed  and  x and  s 
computed,  the  following  specific  acceptance  criteria  for  the  lot  can  be 
applied. 

Accept  the  lot  if 

> L for  i-l , . . . , n (2.2) 

and 


XAn,N-n  5 =*L' 


(2.3) 


Those  criteria  follow  directly  from  the  consideration  that  only  those  lots 
that  have  no  defective  items  in  them  are  acceptable.  The  first  condition 
(2.2)  specifies  that  there  must  be  no  defectives  in  the  sample  and  the 
second  condition  specifies  that  the  sample  must  predict  with  a high  degree 
of  confidence  (90  percent)  that  there  are  no  defectives  in  the  remaining 
portion  of  the  lot. 

2.2  Selection  of  Sample  Site 

In  the  classical  operating-characteristics  up,  •!.  to  indexing  sampling 
plans,  the  acceptance  criteria  and  sample  site  are  chosen  so  that  the  producer 
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. FINITE  LOT  VARIABLES  SAMPLING  PLAN  ACCEPTANCE  FACTORS, 
6 n,  N-n,  FOR  y = 0.90  AND  ACCEPTANCE  QUALITY 
LEVEL  PERCENTILES,  AQLP,  FOR  1 - a = 0.95 


e 

(0.90) 

AQLP 

10.95) 

Lot 

Size 

Sample 

Size 

0 

(0.90) 

AQI.1* 

(0.95) 

Lot 

Size 

Sample 

Size 

e 

(0.90) 

“.594 

14.949 

15 

2 

13.320 

26.143 

25 

2 

15.353 

3.066 

5.460 

3 

5.491 

9.590 

3 

6.256 

1.831 

3.156 

4 

4.104 

6. '32 

4 

4.682 

5 

3.517 

5.534 

5 

4.035 

8.?30 

17.166 

6 

3.174 

4.844 

6 

3.675 

3.61S 

6.38? 

2.935 

4 . 374 

? 

3.438 

2.484 

4.183 

8 

2.745 

4.014 

8 

3.267 

1.680 

9 

2.580 

3.711 

9 

3.134 

10 

2.422 

3.432 

10 

3.02S 

9.619 

IS. 905 

11 

2.257 

3.1SS 

11 

2.931 

4.009 

7.058 

A 'SO 

16 

13.592 

26.676 

12 
1 s 

2.84? 
i 7?n 

3.946 

5.9'0 

3. “01 

5.  160 

3.52“ 

S . 098 

3.395 

1.824 

3.290 

4.606 

3.203 

4.427 

3.129 

4.2" 

3.064 

4.  146 

. 005 

4 . 029 

2.952 

3.925 

2.902 

3.829 

2.855 

3.  Ml 

2.810 

3.657 

2 . ’66 

3.5-8 

.’."22 

3.507 

\M8 

3.42“ 

2.  Ml 

3.353 

2 . 589 

3.2*9 

2.541 

3.  204 

and  consumer's  risks  satisfy  certain  constraints.  Specifically,  the 
producer's  risk,  the  probability  of  rejecting  lots  at  the  acceptable 
quality  level  (AQL)  must  be  less  than  or  equal  to  a and  the  consumer's 
risk,  the  probability  of  accepting  lots  at  the  rejectable  quality  level 
(RQL)  must  be  less  than  or  equal  to  1-y  . For  each  sample  size,  the 
acceptance  region  given  in  (2.1)  is  designed  to  give  the  consumer  lots  of 
high  quality.  It  is  granted  that  (2.2)  and  (2.3)  do  not  lead  to  the 
classical  protection  criteria  at  the  RQL  mentioned  above,  but  they  do  lend 
themselves  to  a consumer1 s- risk  interpretation.  Specifically,  only  those 
lots  that  the  user  predicts  not  to  be  of  rejectable  quality  (one  or  more 
defects)  are  accepted.  Here,  the  consumer  is  "risking"  a faulty  prediction 
at  the  100y  th  confidence  level. 

Taking  this  as  a proper  consumer's  risk  control,  we  are  free  to  choose  the 
sample  site  to  control  the  producer's  risk.  Here  we  cun  closely  Follow  the 
classical  interpretation,  in  that  wo  require  that  lots  with  :ero  defects 
will  be  accepted  at  least  100  (1-a)  percent  of  the  time.  In  order  to  compute 
this  probability, however,  we  must  specify  the  quality  of  the  process  the 
lot  comes  from  as  measured  by  p,  the  average  fraction  defective.  Letting 
<•(•)  lie  the  cumulative  standard  normal  distribution,  n and  u the  mean  and 
standard  deviation  of  the  manufacturing  process,  we  have  that  p="V  ((I.-ol/o). 
If  we  let  k be  the  number  of  defects  in  the  lot,  we  have  that  the  sample 
size,  a,  must  he  such  that: 

I*  (accepting  lot  1 k=0,  p i 1-a  (2rl) 

If  the  number  of  defects  in  the  lot  were  not  specified  as  a condition  in 
(2.-I)  and  only  (2.3)  were  used  as  an  acceptance  criterion,  then  the 
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probability  of  accepting  the  lot  could  be  computed  directly  as  a function 
of  p using  the  non  central  t-distribution.  However,  this  is  not  the  case. 
Furthermore,  an  exact  expression  for  the  left  hand  side  of  (2.4)  is  not 
available  in  general.  Fortunately,  upper  and  lower  bounds  for  the  left 
hand  side  of (2.4)  are  derivable.  These  bounds  prove  to  be  quite  tight  for 
all  values  of  n and  N of  practical  importance  in  the  sampling-plan  applica- 
tion. A theorem  giving  expressions  for  the  bounds  is  proved  in  Section  4 
for  general  k.  If  we  let  £ be  the  lower  bound  and  require  that  n be  such 
that 

P > 1-a  , (2.S) 

then  (2.4)  is  satisfied.  This  is  exactly  what  has  been  done  to  generate 
the  values  in  Table  1.  Specifically,  Table  1 gives  to  each  N and  n tabula- 
ted, the  acceptance  parameter  0^  n N_n  for  y=0.90  and  the  acceptable  quality 
level  percentile  (AQLP)  defined  as 

*p  * (L-lO/0  (2.6) 

with  p = '•'Up)  for  l-<*  = 0.95. 

It  is  interesting  to  note  that  the  upper  bound  I’  for  the  probability  of 
acceptance  is,  for  many  combinations  of  N and  n shown  in  Table  1,  equal  to 
P.  It  is  for  these  combinations  of  values  of  n and  M that  (2.4)  is  an  exact 
equality.  This  occurs  for  all  n less  than  N/2  up  to  N=25,  for  all  n less 
than  N/3  up  to  N=S0,  for  ail  n less  than  N/4  up  to  N=80  and  for  all  n S 22 
for  N=90  and  100.  In  all  cases  the  difference  is  less  than  11  between  P 
and  P for  n S N/2,  and  less  than  0.081  for  n ? N/3.  Values  for  1-0=0.90 
as  well  as  all  values  of  P are  given  in  Fertig  and  Mann  (5). 

Before  wo  proceed  with  a proof  of  the  theorem  referred  to  abo,-e,  we  give 
an  example  of  use  of  the  tables. 
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3.  EXAMPLE  OF  USE  OF  THE  TABLES 


Consider  a lot  of  size  SO.  Suppose  one  wishes  to  reject  the  lot  even  if 
it  has  only  one  defect.  The  binomial  plan,  which  ignores  the  finite  size  of 
the  lot,  would  require  a sample  of  size  114  in  order  to  reject  a lot  ninety 
percent  of  the  time  that  is  2^=100  Cl/50)  % defective.  The  hypergeometric 
plan  would  require  a sample  of  size  45  to  insure  a 90%  rejection  rate. 

Since  it  is  required  to  have  such  high  quality  in  the  accepted  lot,  one 
might  wish  to  use  a plan  based  on  the  assumption  that  the  manufacturing 
process  producing  the  lot  is  very  good.  If  variables  data  are  being 
recorded,  and  it  can  be  assumed  that  these  data  are  normally  distributed, 
then  the  procedure  presented  herein  will  apply.  Specifically,  referring 
to  Table  1,  we  find  that  a sample  size  of  9,  for  example,  will  guarantee  a 
9S%  acceptance  rate  if  the  manufacturing  process  mean  is  at  least  S.147  o 
above  the  lower  specification  limit.  We  would  accept  the  lot  if  the  corres- 
ponding sample,  x^,...,  Xg,  had  no  defects  and,  from  Table  1,  its  mean  was 
at  least  3.632  sample  standard  deviations  above  L.  From  this  latter  criterion 
we  see  why  the  process  mean  must  be  better  than  5o  above  L in  order  to 
guarantee  a 95%  acceptance  rate. 

At  this  point,  it  might  bi  appropriate  to  discuss  a criticism  of  the  above 
procedure,  namely  expecting  the  normality  assumption  to  hold  at  the  extreme 
tails  of  the  process  distribution  (5  o from  the  mean).  In  fact,  however, 
normality  in  the  tails  is  not  as  necessary  as  one  might  first  think  to 
maintain  accuracy  in  the  calculation  of  the  probability  of  acceptance.  The 
normality  assumption  in  the  tails  is  used  only  in  converting  from  the  AQLP, 

'i  ■ Vo,  to  the  process  fraction  defective,  p.  In  fact,  since  p will  tend 
to  oe  very  small,  it  can  be  off  by  many  oraers  of  magnitude  because  of  a 
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lack  of  normality  in  the  tails  and  still  the  probability  of  acceptance  will 
remain  relatively  unaffected.  This  can  be  seen  by  inspecting  the  binomial 
probabilities  in  (4.4)  and  (4.5).  The  important  parameter  in  the  calculation 
of  the  probability  of  acceptance  is  not  the  fraction  defective  per  se,  but 
rather  the  acceptable  quality  level  percentile  (L-p)/o  which  is  directly 
proportional  to  the  noncentrality  parameter  in  (4.4).  It  is  this  parameter 
which  must  be  controlled  by  the  manufacturer  as  well  as  the  process  average 
defective  in  order  to  guarantee  a high  probability  of  acceptance.  It  is 
for  this  reason  that  the  tables  are  given  in  terms  of  AQLP  rather  than  the 
AQL.  Even  if  the  quality  is  not  as  good  as  the  AQLP,  all  that  occurs  is  a 
higher  rejection  rate,  essentially  what  is  occurring  with  the  hypergeometric 
plan  wherein  practically  the  entire  lot  is  being  screened.  On  the  other 
hand,  if  the  process  quality  is  as  good  as  the  AQLP,  a substantial  savings 
in  sample  size  has  been  effected. 

Finally,  we  should  remark  that  requiring  very  high  manufacturing  quality 
levels  is  not  unreasonable  in  most  cases  in  which  lots  would  be  unacceptable 
if  they  contained  even  a single  defect. 
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4.  THE  BOUNDS  ON  THE  PROBABILITY  OF  ACCEPTANCE  OF  A LOT 


In  this  section  we  prove  the  following  theorem  which  gives  expressions  for 
bounds  on  the  probability  of  acceptance  of  a lot  of  size  N on  the  basis  of 
(2.2)  and  (2.3)  when  we  specify  fraction  defective  p in  the  process,  number 
of  defectives  k in  the  lot  and  confidence  level  Y that  an  accepted  lot  has 
zero  defectives,  he  note  that  this  theorem  may  well  be  of  interest  outside 
the  sampling-plan  context.  It  applies  generally  to  truncated-sampling 
applications  and  stratified-sampling  problems. 


Theorem 

Let  T (t|v,6)  be  the  cumulative  non-central  t-distribution  with  v degrees 
of  freedom  and  noncentrality  parameter  6.  Let  Zp  be  the  lOOp  th  percentile 
of  the  standard  normal  distribution.  Let  6 be  the  variables  acceptance 
factor  defined  by  (2.3)  for  a lot  of  size  N and  sample  of  size  n.  Consider 
lots  that  have  k defects  and  come  from  a process  which  produces  fraction  p 
defective.  The  probability  of  accepting  such  lots  satisfies; 


where 


with 


and 


£ <P  (accepting  lot  |k,p)  < P 

P = min  (1,P*) 

P = P*-D, 


P*  = H(0;k,n,N) 

B(0;p,n) 


l-T(v'n~0|n-l,-vnZp) 


(4.1) 

(4.2) 

(4.3) 

(4.4) 


„ . H(0;k,n,N) 

B(0;p,n)  4=1 


ll(H;k,n,N)  = 


1-T(v^n  Og I n-L-l.-ynU  tp)  R(L;p,n)  (4 

(DC) 

( n ) ’ 


p'j  mi-, (4. s) 

B(0;  p,n-{.) 


(4.6) 
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B(£;p,n) 

“ (;)p‘  ci-pj 

n-H 

(4.7) 

g*  =L 

(n-t-l)/  (n-S.)  \ 

1/2 

! 

(4.8) 

\ 

(n-H) (n-l)-nfG2/ 

m = max 

jail  intergers 

< ."(n-1),  1 

(4.9) 

1 

n-l+nO*"  1 

The  values  of  0 in  (4.4)  and  of  0^  in  (4.5),  as  can  be  seen  from  (4.1), 
depend  upon  the  specified  value  of  the  confidence  level  Y as  well  as  n and  N-n. 
Also,  as  will  be  seen  subsequently,  summation  on  S.  in  (4.5)  should  be  taken 
as  zero  if  m=0.  For  this  case  P=P  and, therefore,  the  probability  of  accepting 
the  lot  can  be  computed  exactly. 

To  prove  (4.1),  we  first  note  that  since  the  events  £=0  and  £/0  are  disjoint 
all-inclusive  events,  then 

P {accepting  lot  |k,p)  (4.10) 

= P (accepting  lot  |fc=0,k,p)  P {£=o|k,p} 

+ P (accepting  lot  1^/0, k,p)  P (£/o|k,p) 


Since  we  reject  the  lot  if  i-  SI,  the  second  term  on  the  right  in  (4.10)  is 
zero.  When  k is  given,  the  probability  distribution  of  i is  hypergeometric. 
Therefore, 


P {1=0  | k , p } 


= H(0;  k,n,N) 

When  £=0,  the  decision  whether  or  not  to  accept  the  lot  rests  on  the  value  of 
the  random  variable  L*  = )(-0S.  Thus,  we  see  that 


(4.11) 


P {accepting  lot  |k,p) 

= P {L*>L( 2=0,k,pl  H(0;k,n,N) 


Finding  the  bounds  in  (4.1)  now  reduces  to  finding  bounds  for  the  probability 
that  L*  is  greater  than  L,  given  that  it  is  based  on  a sample  with  no  defects 
coming  from  a lot  with  k defects  produced  by  a process  with  p fraction 
defective.  Since  we  assume  the  items  in  the  lot  are  stochastically  inde- 
pendent observations  from  the  process  and  since  we  are  considering  a series 
of  lots  all  from  the  same  process,  sampling  from  the  lot  and  requiring  the 
sample  to  have  zero  defects  is  the  same  as  sampling  from  the  process  and 
still  requiring  the  sample  to  have  zero  defects.  That  is 

P {L*>L|2=0,  k,  p)  = P {L*>L| 2=0,  p) 


where  the  expression  on  the  right  can  be  computed  assuming  a sample  of  size 
n is  taken  from  a process  of  quality  P with  the  restriction  that  no  defects 
are  in  the  sample.  But  since  the  events  2=2  for  2=0,..,  n are  disjoint 
all-inclusive  event  for  this  procedure,  we  see  that 


P {L*>L|«.=0,p)  - t4'12) 

n P{L*>L|&=2,p)  P(2=2|p} 

-V  

2=1  P (2=0|p) 

In  order  to  proceed,  we  need  the  following  three  lemmas. 


lemma  1 Consider  a sample  of  size  n^  from  a process  of  quality  p.  let 

L*  = Xj-QSj  where  and  Sj  are  the  observable  sample  mean  and 

2 


(Xl=f  V'V  Si  /(nru)) 

P {L*>L |p}  = 1-T  (/TTj  0[n J - 1 , /TTj  z ) 


standard  deviation 
Then 
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Lemma 


2 P {L*>L|p}  < P {L*>L|8.,=0,  P} 

where  is  the  number  of  defects  in  the  sample  of  size 

Lemma  3 P IL*>L|  ,p}  < P {L*"'  >l|x^, . . , x^=L,  ftj  =0,p) 

where  is  formed  by  replacing  all  observations  in  the  sample 
that  are  less  than  L with  the  value  L itself.  (Thus,  , the 
number  of  defects  remaining,  is  zero.) 

* 

The  latter  two  lesimas  follow  easily  if  one  notes  that  is  an  increasing 
function  of  each  of  the  sample  observations.  Thus  restricting  these  obser- 
vations to  be  above  certain  limits  increases  the  probability  of  observing  ar. 
* 

Lj  greater  than  L. 


The  first  lemma  follows  from  the  definition  of  a non-central  t-variate, 
Specifical ly, 


F {L*>L |p)  = P (Xj-eSjHlp) 


4, 


X,-u 


- ^ 


I a 


Sj/O 


> /iTjO 


Since  ==  , yFT  (X-l0/(7  is  N(0,1),  and  s ^/o  is  the  square  root  of  an 

independent  chi  square  variate  divided  by  its  degrees  of  freedom  (iij-1),  the 

j 

lemma  follows  directly  by  recalling  that  t = (z*6 )/x^/\)  is  a non-central  t 


variate  if  z is  standard  normal  and  is  an  independent  chi  square  with 
degrees  of  freedom  v ■ 


Since  the  distribution  of  >’ , given  the  sample  comes  from  a process  of  quality 
p,  is  binomial  with  parameters  p and  n,  we  scr  from  Lemma  1 that  14.12) 
Implies 


(4.13) 


P {L*>LU=0,p}<  (l-T(/n0[n-l,-fntp))/B(O;n,p) 

From  (4.11)  and  (4.13)  we  see  we  can  bound  the  probability  of  accepting  the 

* 

lot  from  above  by  the  expression  for  P given  in  (4.4).  This  establishes 
(4.2). 

In  order  to  establish  a lower  bound  or  the  probability  of  accepting  the  lot, 
one  could  use  Lemma  2 directly  in  (4.11).  A sharper  bound,  however,  will  in 
general  be  given  by  considering  the  summation  in  (4.12).  Let  L*  be  based  on 
an  observable  sample  with  £ defectives,  <L  and  n-£  nondefectives 

Xf  Xn2L'  Let  *2  ■ £t+1V  <"-*>*  s2  * (rikl) 

We  now  seek  to  compute  P(L* '>L|Xj, . . . ,Xg=L,£'=0,p)  defined  in  Lemma  3.  After 
some  algebra,  we  find  that  using  Lemma  3 gives 
P {X-9S>L|*=4,p} 

5 p {(ir  - Itt ) <VL>2  ? ttt-  n^r  ®2  s2 1 V°'p[ 

2 

This  probability  can  be  non-tore  only  if  ^ ^-y>0.  Thus,  H must  be 

strictly  less  than  n(n-l )/ (n-l+rO2) . This  is  the  value  given  for  m in  (4.9). 
Since  the  number  of  defects,  t ,,  in  the  portion  of  the  sample  to  be  used  to 
compute  X,  is  tcro,  we  have  X-,5L.  Therefore,  we  see  that  for  t<«i, 

P {X-t)S>l,|jt=4,p} 

< P j(X2-L)/S2  > 


n0 


1/2 


(n-H)  (n-l)-nW 


lJ2  = °.P  {• 

Since  X,  and  are  based  on  an  observable  sample  of  site  n,  = n-L  that 

has  toro  defects,  we  can  apply  Lemma  2 and  a relation  similar  to  the  first 
half  of  this  proof  (in  which  tho  expression  for  P*  is  derived)  to  find  that 
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P {X-es>L|£=£.,p}  < 


P{X2-0*S2>L|p} 

P{X2=0|p} 


Here,  the  probability  in  the  numerator  on  the  right  is  now 
restriction  on  the  number  of  defects  in  the  sample  of  size 
we  can  apply  Lemma  1 to  get  for  each  £*1,,,.,  m 

(l-T(PrnoJ|n-£-l,-  /frJz  1) 

P {X-OS>L|JUH,p}  S — — E— 

B(0;p,n-£) 


and 


P {X-0S>L|&  £,p)  = 0 


for  £=  m+1 , , , , n. 


From  (4.15)  and  the  definition  of  m,  we  see  that  the  sum  in 
taken  as  zero  if  m=0.  For  m51,  we  see  from  14.11),  (4.12) 


computed  with  no 
n2=n-£.  Therefore, 

(4.14) 


(4.15) 


(4.14)  should  be 

(4.14) ,  and  (4,15), 


that  D is  as  given  in  (4.5).  This  completes  the  proof  of  the  theorem. 


5.  CALCULATION  OF  TABLE  ENTRIES 

The  distribution  of  the  random  variate  R ,,  = (X-Y,  .,  ) / S,  given  by 

n,N-n  1 ,N-n'  & ’ 

(2.1), which  provides  one  of  the  criteria  for  acceptance  or  rejection  of  a 

lot  of  size  N on  the  basis  of  a sample  of  size  n,  is  considered  by  Fertig 

and  Mann  [6] . They  show  that 

P(R  <0)  = 
n,N-n 

— f ll-fOO]^"'1  exp(-z2/2)T  (Ai0|n-1,-  /n  z)dz,  (4.1) 

/2rT  ^ 

where  T( * | V , 6)  indicates  the  cumulative  noncentral  t-distribution  function 
with  degrees  of  freedom  and  noncentrality  parameter  6 and  4>(")  indicates 
the  cumulative  standard  normal  distribution  function  with  mean  zero  and 
variance  one. 

For  each  combination  of  N and  n given  in  Table  1 I’  (R  <0)  = .90  was  solved 

n,N-n 

iteratively  for  0 by  use  of  Muller's  (10)  method,  u multipoint  iteration 
scheme  which  determines  a next  guess  by  using  inverse  parabolic  interpola- 
tion based  on  the  previous  two  iterates  and  their  midpoint.  A first  guess 
for  the  value  of  0,  was  provided  by  the  approximation  suggested  in  Mann, 
Schafer  and  Singpurwal la  (9)  and  investigated  by  Fertig  and  Mann  J6|.  Sixty- 
four  point  llcrmito  Gauss  quadrature  was  used  to  evaluate  the  integral  given 
by  equation  (4.1),  the  function  ♦(•)  having  been  determined  from 
erf  (x)  » 24>  (x/J)-l,  x JO,  with  the  error  function  erf  (•)  computed  using 
a double  precision  IBM  subroutine  and  the  noncontral  t-distribution  function 
computod  using  an  algorithm  of  11. E.  Cooper  (I).  Where  the  combination  of 
values  of  n and  N-n  agree  with  those  of  llahn  (7),  computed  using  the  multi- 
variate t-distribution,  our  computed  values  agree  with  his  to  the  number  of 
significant  figures  given. 
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The  acceptable  quality  level  percentile  AQLP  defined  as  z^  = -(U-L)/cf  , 
with  p equal  to  AQL,  the  acceptable  quality  level,  and  related  to  z^  through 
P “ ^(zp)»  was  determined  iteratively  by  requiring  equality  in  (2.5)  with  £ 
given  by  (4.3).  Muller’s  [10]  method,  discussed  above,  was  used  in  this 
iterative  determination. 
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